Lattice Dynamics and Electron-Phonon Interaction in (3,3) Carbon Nanotubes 
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We present a detailed study of the lattice dynamics and electron-phonon coupling for a (3,3) 
carbon nanotube which belongs to the class of small diameter based nanotubes which have recently 
been claimed to be superconducting. We treat the electronic and phononic degrees of freedom com- 
pletely by modern ab-initio methods without involving approximations beyond the local density 
approximation. Using density functional perturbation theory we find a mean-field Peierls transition 
temperature of «240 K which is an order of magnitude larger than the calculated superconducting 
transition temperature. Thus in (3,3) tubes the Peierls transition might compete with supercon- 
ductivity. The Peierls instability is related to the special 2ki? nesting feature of the Fermi surface. 
Due to the special topology of the (n,n) tubes also a q=0 coupling between the two bands crossing 
the Fermi energy at Icf is possible which leads to a phonon softening at the F point. 

PACS numbers: 63.22.+m, 63.20.Kr, 63.20.Dj, 71. 15. Mb 



During the past 10 years carbon nanotubes have gained 
a lot of attention This is due to their potential for 
applications (e.g. for molecular electronics) as well as 
due to the fact that they allow to study electronic sys- 
tems in one dimension. Thus they offer the opportunity 
to investigate effects like Peierls transition, superconduc- 
tivity, electron-phonon interaction and the interplay be- 
tween them in a low-dimensional system. Despite the 
great interest in these materials progress has been hin- 
dered until recently due to the difficulty to produce car- 
bon nanotubes with well defined radii as well as due to 
the difficulty to determine in detail the structure of nan- 
otubes present in a given sample. Recently, however, 
very promising progress in identifying the structure of 
nanotubes has been made by combining Raman and pho- 
toluminescence measurements 2] . Size selection has been 
also achieved in certain cases. For example, growing nan- 
otubes in zeolite crystals has made it possible to produce 
tubes with a very narrow radii distribution 3] thus allow- 
ing for a detailed comparison with modern density func- 
tional theory (DFT) based calculations 0, 5, 6] . Recently 
superconductivity has been reported for nanotubes with 
radii of 4 A • This immediately rises the question as to 
the origin of superconductivity, the importance of Peierls 
distortions and of electron-electron correlations. Dealing 
with electron-lattice and strong electron-electron inter- 
action in a materials specific way from ab-initio meth- 
ods has not been possible so far even for much simpler 
systems than nanotubes. Thus in the past either the 
electron-electron aspect has been emphasized how- 
ever, restricted to model studies or the electron-lattice 
interaction has been treated in greater detail dealing with 
correlations only on the level of local density approxima- 
tion (LDA) based functionals P.ITo|. 

In the latter case, however, modern DFT based meth- 
ods allow for the parameter free microscopic calculation 
of phonon modes. So far, for nanotubes this scheme 
has been used nearly exclusively for the determination 
of Raman active modes. The results are generally in 



good agreement with available Raman data thus em- 
phasizing the reliability of this approach. Calculation 
of the full phonon dispersion for a given nanotube is 
a much more elaborate task. Based on supercell cal- 
culations the phonon dispersion for a small number of 
tubes could be determined recently [|| [ljj. This ap- 
proach, however, has a certain disadvantage if one is in- 
terested in phonon anomalies and electron-phonon cou- 
pling. Anomalies show up in most cases at phonon wave 
vectors which are not commensurable with the under- 
lying lattice, thus even for approximate treatments huge 
supercells would be needed to study for example a Peierls 
transition by this method. Only in very favorable cases, 
these anomalies appear at high-symmetry points and can 
be dealt with accura tely by the supercell approach (as, 
e. g., for graphite [iflllj)- Also the calculation of super- 
conducting properties in the framework of the Eliashbcrg 
theory requires a detailed knowledge of the phonon dis- 
persion over the whole Brillouin zone (BZ). Especially 
systems which show phonon anomalies require usually a 
very dense mesh of phonon wave vectors which can not 
be obtained with supercell methods. The calculations 
get even more demanding if the phonon anomalies are 
related to special nesting features of the Fermi surface. 
Since in one-dimensional systems the Fermi surface con- 
sists only of isolated points the most extreme case for 
nesting properties is reached here. This requires also a 
very dense k-point mesh for calculating the electronic 
bandstructure and wave functions. Already for simple 
systems like graphene and graphite this effect can be seen 
easily. It is reflected in the high sensitivity of certain 
phonon modes at the K point to sampling effects of the 
Fermi surface as reported recently [lfj, [llj . 

An alternative method is offered by using density func- 
tional perturbation theory (DFPT) which allows for the 
calculation of phonon modes at arbitrary wave vectors 
without relying on large supercells |12(. Using this ap- 
proach we have studied in detail the complete lattice dy- 
namics and electron-phonon interaction for a (3,3) nan- 
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otube. This belongs to the class of small radii tubes 
which have been reported to be superconducting. Our 
study will try to answer the question whether or not 
electron-phonon coupling can explain superconductiv- 
ity in this tube. A special complication arises from 
the Peierls instability which has to show up in all one- 
dimensional metallic systems. However, this has not re- 
ceived much attention in the past since general believe 
was that the mean-field Peierls transition temperature for 
nanotubes is always very low so that it is of no practical 
consequence. These arguments are based on information 
from graphite using the folding concept However, 
this approach breaks down for tubes with small diameter 
which is seen e.g. in the prediction of the wrong ground 
state for small-radii nanotubes 0, 0] . Recently an ap- 
proximate treatment of the lattice dynamics for (5,0), 
(6,0) and (5,5) tubes has been presented which showed for 
(5,0) tube indeed that the estimated transition tempera- 
ture is of the order of 160 K and thus not negligible [l4| . 
However this approach was based on a non-selfconsistent 
tight-binding scheme and did use only an approximate 
treatment of the polarization eigenvectors thus results 
might be questioned. 

In contrast, our DFT-based method allows for a consis- 
tent calculation of electronic states, phonon modes and 
electron-phonon coupling without introducing approxi- 
mations beyond the LDA level. In this paper we present 
fully ab-initio results for phonon dispersion and eigen- 
vectors as well as for the electron-phonon coupling for 
the (3,3) nanotubes using a well tested norm-conserving 
pseudopotential of Hamann-Schluter-Chiang type [lfij . 
We use DFPT in the mixed-basis pseudopotential for- 
malism |16) which has been successfully applied to study 
electron-phonon mediated superconductivity 0,0]. As 
basis functions localized 2s and 2p functions are used to- 
gether with plane waves up to an energy cut-off of 20 
Ry. For integration over the BZ a Gaussian broaden- 
ing scheme is employed. As test of the reliability of our 
phonon approach we have calculated the phonon disper- 
sion for graphene. Comparison with results published 
recently show excellent agreement among the different 
calculations except for the highest mode at the K point 
which is very sensitive to k-point mesh and broadening 
as already mentioned 0, 0] . With a very dense k-point 
mesh of 5184 points in the BZ we found still a fluctuation 
of 1.5 meV for the highest mode when going from a broad- 
ening of 0.05 to 0.2 cV. The authors of Ref. OH tried to 
avoid the complications due to very dense k-point sam- 
pling by increasing the broadening, however, for study- 
ing instabilities due to Fermi surface nesting that is not 
a practical way since any instability will be broadened 
substantially. Thus one of the complication for the cal- 
culations for the nanotube is the requirement of a very 
dense k-point mesh and a small broadening. 

Our calculations for the (3,3) nanotubes were done in 
a supercell geometry so that all tubes are aligned on a 
hexagonal array with a closest distance between adjacent 
tubes of 10 A. The tube-tube interactions are very small 
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FIG. 1: Calculated bandstructure of the (3,3) nanotube. 
Compared are results obtained (a) from the calculated 
graphene bandstructure using the folding method, and (b) 
for the true nanotube geometry. 

[4j. We used 129 k points in the irreducible part of the 
BZ and a broadening of 0.2 eV. The structure was fully 
relaxed and the optimalgeometry agreed very well with 
those given in Refs. 0, |5|- The phonon calculation was 
carried out with the DFPT method. For calculation of 
the electron-phonon matrix elements we even used up 
to 1025 k points. In Fig. ^ we have plotted the band- 
structure close to the Fermi energy and for comparison 
also the bandstructure as obtained by using the folding 
method. Our results agree well with those obtained pre- 
viously 0,0. At k=k F =0.284 ( in units of 27r/a where 
a is the lattice constant of the graphene honeycomb lat- 
tice) we see that two bands are crossing ep- This special 
feature which holds for all (n,n) tubes is important for 
phonon anomalies seen at the T point as will be em- 
phasized later. The dispersion differs most notably from 
those obtained by using the folding technique by a shift 
of the kp value and a change in slope of the two bands 
crossing at cf- This has of course drastic consequences 
for the phonon modes. 

Phonon results are shown in Fig. [21 Again we have 
plotted folding results together with the ab-initio dis- 
persion curves. These ab-initio results were obtained 
by Fourier interpolation on a 16-point grid. The elec- 
tronic wave functions were obtained from a calculation 
with a broadening width of 0.2 eV. In general, the ab- 
initio phonon frequencies are softer than the folding re- 
sults. Furthermore, near q=2k^ one sees phonon anoma- 
lies in certain phonon branches as well as a certain soft- 
ening at the r point which is not present in the disper- 
sion obtained by folding. The T-point softening has been 
seen and studied already in Ref. [3, EH for the G band, 
however, this effect is also present in another mode at 
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FIG. 2: Calculated phonon dispersion curves of the (3,3) 
nanotube as obtained (a) by the folding method from the 
graphene phonon dispersions, and (b) for the true nanotube 
geometry. The latter corresponds to a high effective temper- 
ature of 1096 K. The arrow indicates 2kf (folded back to the 
first Brillouin zone). 



«90 meV with the same symmetry. It is related to the 
fact that these T-point phonons couple the two electronic 
bands crossing right at kp (Fig. [TJ, leading to a nesting 
vector q=0. Therefore, this softening also indicates an 
underlying Peierls instability, however, the correspond- 
ing Peierls distortion does not change the lattice period- 
icity. Note that in our DFPT approach, the calculations 
of phonons include all screening effects due to lattice dis- 
tortions in contrast to other treatments like e.g. tight- 
binding methods (see Ref. [T^jl. 

To study the anomalies near q=0 and q=2kp in more 
detail we have increased the number of k points and re- 
duced the broadening from 0.2 to 0.025 eV. This can be 
interpreted as a variation of the electronic temperature 
from 1096 to 137 K. With finer sampling and reduced 
temperature one of the modes with the anomalous be- 
havior at q=2kp gets unstable. Results for T e ; = 137K 
are shown in Fig. [3J Here only those modes which are 
sensitive to the temperature variation are shown. The 
modes which show anomalous behavior at the T point 
stay stable for all temperatures studied but should even- 
tually go soft at low enough temperatures. Using the 
definition of the Peierls transition temperature (Tp) with 
^2k F (Tp)=0 leads to a lower limit of 137 K for Tp which 
is of similar magnitude as the value given for a (5,0) tube 
in Ref. [l4|. It is important to note that, contrary to early 
estimates of Tp for nanotubes of the order of IK [19j . 
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FIG. 3: Phonon dispersion curves for the two symme- 
try classes which are affected by electron-phonon coupling. 
Shown are results obtained on a fine q grid and for a small 
effective temperature of 137 K. 

this temperature is of sizeable magnitude and thus the 
Peierls instability might compete with the superconduct- 
ing transition. 

To investigate the possibility of phonon-induced su- 
perconductivity, we have calculated the microscopic 
electron-phonon coupling parameters, which are central 
to the Eliashberg theory of strong-coupling superconduc- 
tivity. The coupling constant for a phonon mode qA is 
given by 

A " A = hN(e F )uj A ^ \9t+ q v';kv\ 25 (tkv)5{<ik+ q w), (1) 

" kvv' 

where N(ep) is the density of states at the Fermi energy 
and all energies are measured with respect to ep. The 
electron-phonon coupling matrix element is given by 

9k+qu';ku ~< $k+qv>\5V^ f | (j)kv > 

with SV^ff being the change of the effective crystal poten- 
tial due to a phonon qA and | <f>k v > being the electronic 
wave functions . Since this quantity is very sensitive 
to k-point sampling we have used a grid of 1025 points. 
Due to the selfconsistcnt determination of SV, the ma- 
trix elements include all screening effects in contrast to 
tight-binding approaches. Because the Fermi surface con- 
sists only of the points k=±kp, contributions to the total 
electron-phonon coupling constant A = 1 / N q ^2 qX \\ are 
restricted to q=0 and q=2kp. Roughly 20% comes from 
q=0 while 80% is contributed from 2kp. An estimate 
for T c can be obtained by using the Allen-Dynes formula 
poj . Neglecting the effect of the electron-electron inter- 
action on the pairing, it gives an upper limit of T c =0.833 
wi n exp{— 1.04(1 + 1/A)}, where the effective phonon fre- 
quency is defined as 

u) ln = exp{^^^ X q x ln(u) qX )}. 

q <!>■ 

Using our results for T e ; = 1096K gives A w 0.25 and 
LO\ n w 60.4 meV and as an upper limit T c = 3 K. One has, 
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however, to consider the possibility that the presence of 
the Peierls instability significantly alters this estimate. 
This is because in the expression £Q) for the coupling 
constant \ q \ the renormalized phonon frequencies enter. 
The softening of the critical mode thus produces a strong 
temperature dependence of A, which formally diverges at 
Tp. At the same time u>i n goes to zero. Simulation of 
this mode softening effect using the Allen-Dynes formula 
shows that T c can be enhanced due to the softmode, but 
the obtained maximal value of « 30 K is still significantly 
smaller that Tp. Taking into account a finite electron- 
electron interaction further reduces this value. 

To get a more precise estimate of the Peierls transi- 
tion temperature we have made connection with simple 
model studies which have been carried out in the past 
considering a Frohlich Hamiltonian in random-phase ap- 
proximation |2l| . Extending this approach to the present 
2-band case gives for the the critical mode at q=2kp a 
temperature dependence of the frequency as 

{hw q=2kF {T)Y = A\n{T/T P ), 

where 

A = 2Huj bare y {g„) 2 N„{e F ). 

Here, N„(ep) is the partial density of states per spin of 
band uj bare is the unrenormalized phonon frequency, 
and g v is the electron-phonon matrix element for scat- 
tering processes from kp to -kp within the z^-th band 
induced by the critical phonon at q=2kp. Due to sym- 
metry, only intra-band scattering processes contribute. 
Estimates of the parameters from our ab-initio calcula- 
tions for different effective temperatures lead to Tp «240 
K and Aw (39 meV) 2 . Note that the expression for A does 



not depend on the phonon frequency but only on the 
phonon eigenvector due to the fact that (g v ) 2 ~ l/uj bare . 
Our finding that Tp is an order of magnitude larger than 
T c suggests that for the (3,3) nanotube the Peierls insta- 
bility dominates over the superconducting one. 

In summary we have presented here a fully ab-initio 
calculation of the lattice dynamics and the electron- 
phonon coupling for the (3,3) carbon nanotube. With- 
out relaying on simplifying approximations the Peierls 
transition could be seen and a Peierls mean field tran- 
sition temperature Tp « 240 K was predicted. Cal- 
culating the electron-phonon coupling using the same 
scheme resulted in a superconducting transition temper- 
ature which was one order of magnitude smaller than 
Tp. This makes it not very likely that superconductiv- 
ity based on the electron-phonon mechanism is present 
in (3,3) nanotubes. This does, however, not exclude the 
possibility that phonon-mediated superconductivity ex- 
ists in the case of the (5,0) nanotube, which is another 
member of the class of a 4 A tubes for which supercon- 
ductivity has been observed experimentally Q. Effects 
which have been discussed here are not accessible with 
the folding method since they depend sensitively on the 
curvature of the nanotubes, the detailed bandstructure 
and nesting features. Since in our approach arbitrary q 
points are treated on the same footing we could easily 
see that the q=0 anomaly is only a special case of strong 
electron-phonon coupling. In contrast to former studies 
we could show that this effect is not restricted to the G 
band but should be seen in another mode, too, belonging 
to the same symmetry class. 
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